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Abstract 

This is a brief review of recent progress in constructing solutions to the matrix model Virasoro 
equations. These equations are parameterized by a degree n polynomial W„{x), and the general 
. solution is labeled by an arbitrary function of n — 1 coefficients of the polynomial. We also discuss 

I in this general framework a special class of (multi-cut) solutions recently studied in the context of 

■ TV = 1 supcrsymmctric gauge theories. 

(N 
O 

:^ 

^ Introduction. It was realized in the beginning of nineties that matrix models partition functions 
Q-i, typically satisfy an infinite set of Virasoro- like equations [1, 2]. These were nothing but Ward identities 
(— I [ (Schwinger-Dyson equations) which mainly fixed matrix model partition functions (because of the 
' topological nature of matrix models [3], the Ward identities were restrictive enough). Moreover, it 
. turned out that one of the most technically effective ways to deal with matrix models was to solve 
rS [ these Virasoro equations (they are also sometimes called loop equations) [4, 5, 6, 7, 8]. 

' At early times of matrix models one usually dealt with Virasoro equations describing relatively 

simple "phases" so that the equations had unambiguous solutions. An interest to more complicated 
phases of matrix models has revived after G. Bonnet, F.David and B.Eynard [9] proposed to deal with 
the multi-cut (or multi-support) solutions (known for long, [10, 11, 6, 12]) in a new way: releasing 
the tunneling constraint [11]. Their approach was later applied by R.Dijkgraaf and C.Vafa [13] to 
description of low energy superpotentials in TV = 1 SUSY gauge theories, [14]. 

More concretely, the authors of [14] considered the M = 2 SUSY gauge (Seiberg-Witten) theory 
in special points where some BPS states become massless. Therefore, these states can condense in 
the vacuum which breaks half of the supersymmetries (leading to TV = 1 SUSY) and gives rise to a 
non-zero superpotential. Values of this superpotential in minima are related to the prepotential of 
a Seiberg-Witten-like theory. In turn, R.Dijkgraaf and C.Vafa associated [13] the prepotential with 
logarithm of a partition function of the Hermitean one-matrix model in the planar limit of multi-cut 
solutions (it was later proved in [15]). 

In fact, actual definition of the multi-cut partition functions is a separate problem. For instance, 
one could simply define them as (arbitrary) solutions to the corresponding Virasoro equations (D- 
module point of view). Then, one may ask what is special about the concrete Dijkgraaf-Vafa (DV) 
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solutions. In the present paper we make further steps in this direction, discuss the space of all solutions 
to the Virasoro equations in the multi-cut phase and show that the DV solutions form a basis in this 
space. They are distinguished by a special property of isomonodromy that allows one to associate with 
these Seiberg-Witten-like systems a Whitham hierarchy [15, 16], the corresponding partition function 
having a multi-matrix model integral representation [9, 17]. 

Hermitean one-matrix model. Hermitean one matrix model is given by the formal matrix 
integral over N x N Hermitean matrix M 

Zw{t) = -^ f DMe^p-(-TrW{M) + TrY,tkMA (1) 

Here W{x) is an arbitrary function that we usually assume to be a polynomial of degree n + 1, 
Wnix) = Y^^'^^ Ti;x'' and the constant g'^ controls the genus expansion. This integral still needs to 
be defined. One possibility is to substitute it with its saddle point approximations [9, 17]. Different 
saddle-points M = Mq are given by the equation W'{Mo) = 0. If the polynomial 

n 

W'ix) = l[{x-ai) (2) 

i=l 

has roots a^, then, since Mq are matrices defined modulo L''(A^)-conjugations (which allow one to 
diagonalize any matrix and permute its eigenvalues), the different saddle points are represented by 

Mo = diag{ai, . . . ,ai;a2, . . . ,02] ■■■ ■,an,...,an) (3) 

with ai appearing iVj times, iVj = N. In fact, there is no need to keep these Ni non-negative 

integers: in final expressions (like formulas for the multidensities and prepotentials) they can be 
substituted by any complex numbers. Moreover, can depend on Tfc (i.e. on the shape of W{(f))) 
and g. 

Now, using at intermediate stage the eigenvalue representation of matrix integrals, one can rewrite 
[9, 17] the matrix integral (1) over N x N matrix M as n-matrix integral over Ni x Ni matrices 
Mi (each obtained with the shift by af. just changing variables in the matrix integral (1)), which is 
nothing but the DV solution [13] 



z'^'''^\t\Mo) ~ j n^M.exp ^T¥t»MM n 



n 

2N,Ni 



i=l \ i,k j i<j 



X exp 2 > {-f- — 7—1 



(4) 



Trj Mf Tr^- M\ 
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The variables t), ' are given by the relation 

oo/n \n/oo \ 

E *4 E + = E E (5) 

k=0 \i=l / i=l \k=0 / 

with arbitrary Ni x Ni matrices Mj. 

Virasoro constraints. The other possibility is to observe that (1) satisfies the infinite set of 
(Virasoro) equations (=Schwinger-Dyson equations, = Ward identities) [2] 



LmZw{t) = 0) m > — 1 



a,b>Q 



d ^ ^ (6) 
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dZw . dZw n wi, n I 1 f'7\ 

W + ^ = ° Vfc = 0,...,n + 1 (7) 

and call any solution to these equations the matrix model partition function. Then, the partition 
function is not a function but a formal D-module, i.e. the entire collection of power series (in t- 

variables), satisfying a system of consistent linear equations. Solution to the equations docs not 
need to be unique, however, an appropriate analytical continuation in t-variables transforms one 
solution to another, and, on a large enough moduli space (of coupling constants t), the whole entity 
can be considered, at least, formally as a single object: this is what we call the partition function. 
Naively different solutions are interpreted as different branches of the partition function, associated 
with different phases of the theory. Further, solutions to the linear differential equations can be often 
represented as integrals (over spectral varieties), but integration "contours" remain unspecified: they 
can be generic chains with complex coefficients (in the case of integer coefficients this is often described 
in terms of monodromies, but in the case of partition functions the coefficients are not restricted to be 
integer) . A model of partition function is an integral formula which has enough many free parameters 
to represent the generic solution of the differential equations in question. 

A familiar example that could clarify these notions is provided with the cylindric functions. Their 
defining equation is 

X^dl + Xdx + - 2.2)] z^{X) = (8) 

and an integral representation is 

ZJx) = — I e-'^^^^+'''^de (9) 
27r Jc 

The model is given by the generic linear combination of two contours, say, chosen as in 8.423 of [18] 
(this choice fixes as the basis the Hankel functions). 

Loop equations. Another form of the Virasoro equations (the loop equation) is produced by 
rewriting the infinite set of these equations through the unique generating function of all single trace 
correlators 



p^'\z\t) = V{z)J^, V(z) = ^-^^, J^^gHogZw (10) 
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fc>0 ^ 



Introducing notations v{z) for J2k ^kZ^ and [...]_,_ ( [•••]_) for the projector onto non-negative (negative) 
degrees of z, one obtains the loop equation [4] 

W'{z)p^'\z\t) = {p^'Hmf + fw{z\t)+9''V{z)p^'\z\t) + [v'{z)p^'\z\t)]_ (11) 

fw{z\t) = [W'(z)p(i) {z\t)] ^ = R{z)J' (12) 

In order to consider (connected) multi-trace correlators, one needs to introduce higher generating 
functions (also named loop mean, resolvent etc) which we call multi-density 

p(-)(2i,...,z^|t) = V(zi)...V(z„)JP (13) 

In what follows, we consider solutions to the Virasoro equations (6) as a formal series in t-variables, 
as well as a series in the coefficient (genus expansion): 

p>0 
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Main results. 



Here we are going to review briefly the main results of the papers [7] and [8], where we 
defined the matrix model partition function as any solution to the Virasoro equations 
(6). In forthcoming paragraphs we briefly comment on these results. 

• Any solution to the Virasoro constraints (taken as a formal series in f-variables 
and in genus expansion) is unambiguously labeled by an arbitrary function of n of 
n + 2 T-variables: the bare all genera prepotential. 

• There is an evolution operator that generates from the t-independent bare pre- 
potential the matrix model partition function which depends on t-variables and 
satisfies the Virasoro equations. This evolution operator does not depend on the 
choice of the arbitrary function, but only on T- and t-variables. 

• One may invariantly define "the occupation numbers" of [9, 13] as eigenvalues 
of operators constructed from the evolution operator, formula (35) below. The 
corresponding DV solutions are described as eigenfunctions of these operators. 

• (Conjecture 1) The evolution operator can be completely expressed in terms of 
the unique operator 



y ^ ^JW'ixr - Vi?(x), Rix) ^ - ^ (a + 6 + 2)r„+5+2X«— (15) 

a,b=0 

its derivatives and W'{x). 

• (Conjecture 2) The evolution operator is constructed as a formal series in t with 
operator coefficients acting on the bare prepotential. These coefficients generate the 
full matrix model correlators. These operator coefficients are related to operator 
multidensities (13) exactly as the full correlators are related to the connected 
correlators, only an appropriate ordering prescription should be applied. This 
relation is universal, i.e. is the same for the Gaussian (quadratic) and non- 
Gaussian potentials. 

• (Conjecture 3) The ordering used in the previous conjecture is not uniquely 
defined. 



Solving the Virasoro constraints. In order to convert (11) into a solvable set of recurrent 
relations, we expand p^^\z\t) in powers of and t's 



p,m>{) 

fw{z\t) = ^ j> ■■■ ^v{zi)...v{Zm)fw'^^^\z\zi,...,Zm) 



p,m>0 

In this way, we introduce the full set of multidensities p^"^\zi, . . . , Zm) and auxiliary polynomials 
j^\m+i)^^^^^^ ^^^^ (which distinguishes between different phases for a given W{z)) at zero t's. 
Acting on eq.(ll) with the operator V{zi) . . . V{zjn) and putting all tfe = afterwards, we obtain 
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a double-recurrent (in p and m) relation for the multidensities p^\'^^ 

{Z,ZI,..., Z^) f'^"'+^{z\zi,. ..,Zm) = 
q mi+m2=m (1"^) 



1=1 



Together with (12) this relation is enough to find explicit expressions for all the multidensities through 
W{z) and f^^^\z). In fact, the latter polynomials (all of degree n — 1) are not independent, since for 
m = 0, 

f^^'\z)=RF(P^[T] (18) 

expresses all the /-polynomials through a single function of T (i.e. of W{z)) and g, the prepotential 

at i = 0, 

oo 

J^it = 0,g)= F[g, T] = E Q^'F^'^ [T] ^ Z[g, T] (19) 

The operator R here is given in (15) and contains derivatives with respect to the T-variables. We 
call such operators check operators and denote by the "check" sign to distinguish them from the hat 
operators, which contain t-derivatives and are denoted by "hats". 

Note that the Tk dependence of F[g,T] is not arbitrary, since the first two (L_i and Lq) Virasoro 
constraints are linear in derivatives and can be consistently truncated to i = and then allow one 
to express two derivatives, say, dF/dTn+i and dF/dTn through dF/dTi with / = 0, . . . ,n — 1. As a 
corollary, the partition function can be represented as 

ZiT)\^^, = Z[g,T] = I dfez(fc;7?2,...,7?n;?i)e^("'"') (20) 

with an arbitrary function 2; of n arguments {k,r]2, ■ ■ ■ ,r)n) and h. Here the L_i-invariant variables 
are used, 

W = — logr„+i, X = To + . . . ~ Tjn+l, 

n + 1 



_ , . k{k-2y. '^ (n + l)^(n-0! ^k-l-i^i W"^ 
As an immediate corollary of (17), we obtain for p = and m = 



(21) 



p{0|i)(^)^EM_ll(£) (22) 

with 

y\z) = {W\z)f-Af^^\'\z) (23) 

Evolution check operator. The basic claim is that, for any W{z), there is an evolution (check) 
operator [/vi^(t), which converts any function Z[r] of Tq, . . . , T„_i (with prescribed dependence on T„ 
and Tn+i) into Zw{t) = Uw{t)Z[T], which satisfies the Virasoro constraints, Lm{t)Zw{t) = 0, m > 
— 1. This means that the evolution operator is the same for any values of the arbitrary parameters 
/ (or for any function Z[T]) once W{z) is fixed and that "orbits" of the evolution operators are 
completely parameterized by W{z). Moreover, if Z[T] = Y^a'^aZ^"'\T], then Z]y{t) = J2a'^aZy^\t). 
This means that one may arbitrarily choose a basis in the space of all solutions, with the evolution 
not disturbing the expansion of any solution w.r.t. this basis. 
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make 



use 



One may construct the operator Uwit) with the following procedure. For given T's, we majvc 
of the Virasoro constraints LmZ{t) = and their multiple t-derivatives to recurrently express 

with all < Ij < n 



dtk^...dtkp 



t=0 



with all < < oo through g'^^ dT^^^dTi 



dPZ 



dtki ■ ■ ■ dtkp 



t=o 

.ll...lp-\-s 



t=0 



3 (.1 ^ 

/l,...,/p-|-s 



t=0 



(24) 



This is a straightforward procedure, and the sum over s is finite, from to the integer part of 

but with a,b < k — n — 1, further, the expression for 



the expression for ^ 

contams p.. p.. „, 

j_Q ataOtiOtc 



contains ^ 



this 



Let us now define the operators 



ftaCifc (=0 7 - _ -1 

time with a,b,c < k — 2n — 2 and so on. 



D 



,{p) 

/c . . . kp 



E 



(25) 



and construct 



(26) 



0<h,.-;lp+s<n—l 

the evolution operator tJw{t) as a series in these £>-operators 

tJw{t) = 1 + tkb^^^ + \tktiDfi> + \tktitmbf;>^ + ... 

The fact that, for any Z\T], 

L^{t)Zw{t) = Lm{t)Uw{t)Z[T] = (27) 

or, simply, that 

\ k a+b=m / ' \ a+b=m ) 



+^ E I ^1 + ^2 A^^i^,,, + E ^T^^^klmM, + E ^Iti,^. 1 + 



n(2) 



n(3) 



n(4) 



a+b=m 



(28) 

is equivalent to vanishing of all the linear combinations of operators in brackets, and these are the 
characteristic equations for the £>-operators. 



Basis in the space of all solutions. One can now cleverly choose some basis in the space of all 

solutions. Note that the DV solutions do form such a basis. To have them written in a more clever 
way, one may present the contribution of a particular extremum Mq (labeled by the set of Ni) above 
in the Givental-style decomposition form, expressing it through the product of n Gaussian partition 
functions {Z^{t\N) given by the N x N matrix integral with quadratic Wg{x) = Mx^), with its own 
Ni and Mi = W"{ai) = Uj^Mi - 



n n 



Matr) ,^„, _ nr=ie-^'^("')Vol^(^,^) ^ ^2A..A., ^ J{d,J{ Z^\t^^ |iV,) (29) 



Z^^^'^\t\M,) = 



Vol 



[/(AT) 



i<j i=l 



with operators 



' [ k,J ^;'km dt^dt^^J' 



(30) 
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a. 



ai — aj, and 



(31) 



fe>3 ^"fe 

(T4/(^)(x) = d^W{x)). 

These DV solutions have a series of properties that basically have much to do with isomonodromic 
deformations (by W{x)) and are in charge of Whitham integrable systems behind the DV solutions 
[15, 16], see the end of the next paragraph. 



Evolution operator as a function of y. So far we basically considered the "connected" corre- 
lation functions, p^'^^'^zi, . . . ,Zra',g)- The other possibility is to deal with "full" correlation functions. 



if(»(zi, ...,zm;g) = Zwit; gr'V{zi) . . . V{zm)Zw{t; g) = £ g'^^^-"'^ K^^^^^zi, ...,Zm) 



p=0 



(32) 



These are related by 



EE E E 



-,2(P1 



K^-\^\zi,...,Zm;g) = 
( 

E 



■+Pv-l>) 



<T A;=l i/i,...,i/fc=lpi,...,pi/=0 



1 



(pi|mi 



\^a{l)-i ■ ■ ■ 1 ^a{mi))H l^cr(mi+l); • • • ) ^cr{m2)) • • • F V^(T(m-m^+l) i • • • ) ^a{m) ) 



z^i!(mi!)^i ...z/fe!(mfc!; 

\ m=vimi+...+Vkmk 

{pv\rhi,)i 



(33) 



Our task is to express the correlation functions defined in (32) and (13) with the help of hat- 
operators through the action of check-operators. In fact, this task is already solved for (32): one can 
easily see that 

X(-H(zi,...,z„)= 



1 



■■■Zm 



(34) 

K(» izi,...,Zm;g)= Z{T- g)-^k^-\^) (zi, . . . , Zm)Z{T- g) 

Moreover, manifest examples that can be found in [8] show that these quantities are expressed through 
y (15), its derivatives and W'{x). 

Note that one may now invariantly define the quantities that emerged in the DV solutions. Si = 
To this end, one should introduce the "occupation number" operators (in fact, these operators describe 
the monodromy of K^'^^\z))^ 

gSi = Ni= (f k^-\^\z)dz (35) 

Then, Si are nothing but the eigenvalues of this operator. This is analogous to the condition ^ ~ 
which one usually adds to (6) in order to describe matrix integrals. Now one can define the DV 

solutions as eigenfunctions of the set of operators Si. Therefore, Si being eigenvalues, by definition, 
do not depend on T^'s. This exactly expresses the isomonodromic properties of the DV solutions, see, 
e.g., [6, 15, 12, 19]). 



Relations like 



K 



dSi 



Ni, 



dNi 



2 r 

^ 9 Oij 



J0\1) 



should serve as operator counterparts of the celebrated 

Si = 

dFpv 

dSi 

for the particular DV solution. However, one should be careful about regularization, higher-loop corrections etc. 



J0|1) 
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Check p-operators. The connected correlators p are more "fundamental" than the full K. There- 
fore, it is natural to wonder if one can find check-operator analogues of p's, once we see that check- 
operator counterparts of K do exist and can be of some use. This means that we would like to 
put 

oo 

m\z-g) = Y,9''''-''p^\'-\z-g) (36) 

In this way, one gets rid of the terms with explicitly present operators R and prepotentials F'^p\ 
the relevant check-operators p^ l^) are expressed through y{z;g) only (with the single exception of 
pC^W (^z; g), which also contains W'{z).) Thus, the check-operator k('\p^ is a polynomial in W' of 
degree p. Instead, the g dependence is now distributed between explicit factors like g^^~^ and an 
additional ^-dependence of y{z;g). This, however, allows check-operators p^\'^^ to look exactly the 

same (modulo ordering) as the corresponding Gaussian multidensities p^Q^\ which are all expressed 
through yc only. 

In [8] we suggested a hypothesis that eq.(36) is indeed true in all orders in g"^ and, moreover, 
similar expansions hold for all ^(»(zi, . . . , z^;^): they can be all expressed through 
multilinear combinations of check operators p^\'^\ which (for {p\m) ^ (0|1)) depend only 
on y{z; g) and its z-derivatives in exactly the same way as p^Q^^ depends on yciz). However, 
even to formulate this hypothesis, one needs to introduce some ordering prescription for products of 
check-operators which is not, as usual, unique. Different ordering prescriptions lead to different explicit 
formulas for p, and our hypothesis states that there exist orderings, when these expressions contain 
y, it derivative and nothing else, except for a few W' , see [8] for more possible choices. 

In principle, when introducing p-operators, we have different possibilities of definition, preserving 
one or another kind of their relation to ^'s. They could be defined similarly to (34) from (13), so 
that eq.(38) below becomes an equality. However, it appears more interesting instead of preserving 
the equations, to require for p('l^) to be the same (up to ordering) as the Gaussian functions p^}^"^ ■ We 
can construct recursively an operator modification of expression (33). Since the recurrent equations 
for K are linear 

W'(z)i^O+i) (^, . . . , _ i2(^)i^(-l-)(^,, . . . , 

d K(-\^\z, Zi,...,Zi,...,Z,n)- K(-\^\zi, ...,Zm) 2 r.(.\m+2) , . (3'^) 

+ 2.^^ '^^^ =9K^^ + >{z,Z,Zl,...,Zm) 

they coincide with the equations for operators K. The equations for functions p arc not linear. Thus, 
for the operators p we should choose the order in which different p stand in the products. 
Note that, with our definition, 

pg"^^ {zi,...,zm)^ Z{T- g)-'p^^^ (zi, . . . , zm; g)Z{T- g) (38) 
and, in variance with the l.h.s. of (38), its r.h.s. is still ^-dependent. 
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